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i Abstract 

Generalizations of the complex number system underlying the 
mathematical formulation of quantum mechanics have been known for 
i— i some time, but the use of the commutative ring of bicomplex numbers 

! for that purpose is relatively new. This paper provides an analyti- 

J> cal solution of the quantum Coulomb potential problem formulated 

in terms of bicomplex numbers. We define the problem by introduc- 
ing a bicomplex hamiltonian operator and extending the canonical 
commutation relations to a form like [Xi, P&] = i-^^a^I, where £ 
[- — is a bicomplex number. Following Pauli's algebraic method, we find 

the eigenvalues of the bicomplex hamiltonian. These eigenvalues are 
also obtained, along with appropriate eigenfunctions, by solving the 
extension of Schrodinger's time-independent differential equation. Ex- 
• i-h amples of solutions are displayed. There is an orthonormal system of 

^ solutions that belongs to a bicomplex Hilbert space. 
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1 Introduction 



It is generally believed that the best justification of a physical theory, no 
matter what the nature of its mathematical formalism, rests on the agreement 
of its predictions with experiments and the internal consistency or elegance 
of the theory itself pQ. Views on the mathematical formalism range from 
pure instrumentalism to the idea that the formalism itself is real [2]. 
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The mathematical formalism of quantum mechanics has been studied 
thoroughly P, HJ |5]. Not everyone agrees on the set of postulates neces- 
sary to build the foundations of a coherent quantum mechanics [61 [7]. The 
Hilbert space structure of the set of quantum states, however, seems to be 
uncontroversial. 

Hilbert spaces used in quantum mechanics are defined over the nonordered 
field of complex numbers C. Complex numbers make up a division algebra 
richer than their real subset, and are deeply connected with superposition of 
quantum-mechanical amplitudes. Unitary representations of Lie groups, fun- 
damental tools in the quantum theory of symmetry, require complex numbers 
in an essential way [SJ E] • 

Quantum mechanics postulates that the only possible results of the mea- 
surement of a dynamical variable are the eigenvalues of the corresponding 
self-adjoint operator acting in the state space. As eigenvalues of self-adjoint 
operators are real, this is a form of correspondence between the quantum 
and classical world descriptions. 

If complex numbers are so appropriate to describe the quantum world, one 
can ask whether generalizations of that number system might do equally well 
or even better. The noncommutative field of quaternions has already been 
investigated from that point of view [10]. More recently, attention has turned 
towards the commutative ring of bicomplex numbers [TT| [T2j. This paper is 
part of a program of extending the quantum mathematical formalism to that 
algebraic structure, which is neither a division algebra nor an absolute- valued 
algebra over the real numbers [131 H2 EE1 EH [T71 [18] . 

In sect. [2j we review some of the algebraic properties of bicomplex num- 
bers and construct an infinite-dimensional bicomplex Hilbert space made 
up of square-integrable bicomplex functions. Section [3] defines the bicom- 
plex generalization of the quantum-mechanical Coulomb potential problem. 
Eigenvalues of the corresponding hamiltonian are obtained through Pauli's 
algebraic method. Section [| is devoted to obtaining eigenfunctions of the 
bicomplex Coulomb potential hamiltonian in the coordinates basis. To our 
knowledge, this is the first time that this has been done with an algebra larger 
than C. Graphical representations of some functions of interest are shown in 
sect. [5] In sect. [6] we check the consistency of some of the assumptions made 
and show that the Coulomb eigenfunctions live in a bicomplex Hilbert space. 
We conclude in sect. [3 
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2 Bicomplex numbers and functions 



In this section we briefly summarize relevant algebraic properties of bicom- 
plex numbers. More information and proofs can be found in [TTJ [TJ]. We 
define the concept of a square-integrable bicomplex function. Such functions 
are then used to construct an infinite-dimensional bicomplex Hilbert space 
which, it turns out, will be an appropriate arena for bicomplex Coulomb 
potential eigenfunctions. 

2.1 Algebraic structure 

One way to define a bicomplex number a is by writing 

a := o^e-L + a%e 3 , (1) 

where and belong to C(i 1 ), the field of complex numbers. The caret 
notation [15] is used to label complex components of bicomplex numbers, 
thereby avoiding confusion with other kinds of indices. The imaginary bi- 
complex units e 1 and e 2 satisfy the remarkable properties 

e i = e n e 2 = e 2> e i + e 2 = 1 an d e i e 2 = = e 2 e 1 . (2) 

We call {e 1 ,e 2 } the idempotent basis of bicomplex numbers. With the ad- 
dition and multiplication defined in the obvious way, the set of bicomplex 
numbers T forms a commutative ring with unity. Properties (2| greatly sim- 
plify bicomplex algebraic calculations and make definition (IT), among other 
equivalent choices, a very useful one for our purposes. 

If a± = 0, then ^ implies that ae 1 = 0. In fact any bicomplex number 
a for which either = or = is a zero divisor. The set of all zero 
divisors is called the null cone and is denoted by MC. 

Define j, the imaginary hyperbolic unit, as j := e 1 — e 2 . Then j 2 = 1, 
e-,^ = (1 + j) /2 and e 2 = (1 — j) /2. Substituting this in Q, we get the 
hyperbolic representation of a as 

= jo a + a j + jo a -o S j j=:![o+B>j) (3) 

which in turn means that = x a + y a and = x a — y a . If a± and 
are both in IR, we call a a hyperbolic number. The set D of all hyperbolic 
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numbers is obviously a subset of 

There are several ways to define conjugation in T, but the one most useful 
for our purposes is the following. We define afi asafe 1 + (^e 2 , where the 
upper bar denotes the usual complex conjugation. Clearly, 



or a = | a j | e x + \at%\ e 2 , (4) 



where | | is the standard real norm of complex algebra. 
The real norm of a bicomplex number a is defined as 
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M : = ^7fVM 2 + l«2l 2 - ( 5 ) 



If a is hyperbolic, then 



a 



y/x 2 a + y 2 a := yjRe(a) 2 + Uy(a) 2 . 
One can show [H] that for all a, (3 G T and z G C(i 1 ), 

|a| > 0, \za\ — \z\ \a\ , \a + f3\ < \a\ + \f3\ and \a/3\ < y/2 \a\ . 

Since T possesses zero divisors and since the norm of a product of bicomplex 
numbers is not in general equal to the product of their norms, the algebraic 
structure (T,+,-,| |) is neither a division algebra nor an absolute-valued 
algebra over the real numbers. 



2.2 Normed function space 

Defining a bicomplex function f of q bicomplex variables as a g-tuple infinite 
positive-integer power series with bicomplex coefficients, one can show that 

/(m) = /t(mi) e i + h(^2) e 2- ( 6 ) 

The notation f(/-i) means that / depends on q bicomplex variables and 
each fs (s = 1, 2) is a C(i 1 ) complex functions of the q complex variables fiis- 
We say that / in (|6| belongs to the null cone if either or fa is zero. We 
call / a hyperbolic function if /-j and are both real. 

1 Note that (— i x j) 2 = —1, so that — i-J has the properties of an imaginary unit. It is 
usually called i 2 . 
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As a particular case, if all real we simply have 

/(M) = /i(M) e x + e 2 . (7) 

We say that / is a bicomplex square-integrable function if and only if the 
are both square-integrable functions, that is, 

J\fM\ 2 dfi<oo (8) 



for s — 1 and 2. Here d^ is the Lebesgue mesure on M. q |19j . We denote 
by J-g the set of bicomplex square-integrable functions of q real variables. It 
can be shown that with standard addition and multiplication, T q makes up a 
T-module. This module is explicitly denoted as (J- q , T, +, ■) and it obviously 
has infinite dimension. 

For any f,g e JF q , the following binary mapping takes two bicomplex 
square-integrable functions and transforms them into a unique bicomplex 
number: 



(f,g) 



J /t( M ) g fr) d ^ = e s J h(»)gs(») dj*. (9) 



If we identify functions that differ only on a set of measure zero, the binary 
mapping ^ satisfies all the properties of a scalar product. Explicitly, 

1. (f,g + h) = {f,g) + (f,h); 

2- (f,ag) = a(f,g); 

3. {f,g) = {g,h)h 

4. (/, /) = if and only if / = 0. 

The functions / and g are orthogonal if their scalar product vanish. We say 
that / is normalized if (/, /) = 1. 

With @, one can define an induced T-norm on T q as 



Making use of results proved in [17J, it is not difficult to show that the 
structure (J-g, T, +,-,(,), || ||) is a bicomplex Hilbert space. 
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3 The Coulomb potential problem 



In standard quantum mechanics, the hamiltonian associated with the Cou- 
lomb potential is given by 

a = J-p° - z 4, (ii) 

where fi, e 2 and Z are positive real numbers and 



P 2 :=Pf + Pj + P*, R:=sJXl + Xl + Xl (12) 

The six operators Xi and Pk are self-adjoint and satisfy the commutation 
relation [Xi, Pk] = '^hdikl. 

The quantum-mechanical Coulomb problem consists in finding the eigen- 
values and eigenvectors of H, that is, it consists in solving the equation 
H\4>e) — E\iPe). Its most important application is the determination of the 
energy levels and state vectors of a hydrogen atom or hydrogen-like ion in 
its center-of-mass frame. The Coulomb problem is one of the few analyti- 
cally solvable problems of quantum mechanics. The solution can be obtained 
both by an algebraic method that goes back to Pauli [2111 EEJ [22] , and by a 
differential equation method that goes back to Schrodinger [2"3"l |2U [22]- 

3.1 Statement of the bicomplex problem 

The quantum Coulomb problem will now be formulated in terms of bicomplex 
numbers. The crucial step consists in extending the canonical commutation 
relations the way it was done for the harmonic oscillator [15] . We proceed by 
making a set of assumptions, from which we will derive a number of properties 
satisfied by eigenvalues and eigenkets of the hamiltonian. The consistency 
of the assumptions will eventually be checked through the explicit solutions 
obtained. So here are our assumptions: 



Seven bicomplex linear operators X^, Pk and H, related by (11) and 



(12), act in a T-module Ai. Elements of M. are called kets and are 
generically denoted as \ip). Operators and kets can be decomposed in 
the idempotent basis as 

Xi = Xne-L + X i2 e 2 , (13) 
|^> = |^>ie 1 + |V> a e a . (14) 
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Here X is : = e s Xj and \i[)) s := e s |^). We say that Xi belongs to the 
null cone if either Xn or X^ = 0. The same applies to Pk and H. 
Similarly, e AfC if \if)i or \if)2 = |0). The bold index notation 
is such that any quantity affected by a bold subscript s is unchanged 
when multiplied by e s . 

b) The operators Xi and Pk are self-adjoint with respect to a bicomplex 
scalar product to be specified explicitly. The scalar product has to 
satisfy properties analogous to the ones enumerated after eq. Self- 
adjointness is denoted as Xi = X* and Pk = PI- 

c) The scalar product of a ket with itself belongs to D + := {a^e 1 + a$e 2 '■ 
Qfj, 02 > 0}. 

d) [Xi, Pk] = i-Jldik^I-, where £ G T is not in the null cone, / is the identity 
operator on M, ti is Planck's reduced constant and 5ik is Kronecker's 
delta. 

e) There are eigenkets \iPe) of H which are not in the null cone and whose 
corresponding eigenvalues E are not in the null cone. 

f) Eigenkets \4>e) corresponding to a given eigenvalue E span a finite- 
dimensional T-module. 

g) Two eigenkets \ipEi), \^e) of H, not in AfC, with (Ei — Ej) not in AfC, 
are orthogonal. 

Assumption (a) introduces the bicomplex generalization of the position, 
momentum and energy operators. With (b) we impose, as in the standard 
case, the self-adjointness of Xi and Pk- The third general property of the 
scalar product, stated after eq. ([9]), implies that (|?/>), \?p)) E D. The more 
restrictive assumption (c) is added so that a ket \ip) ^ AfC can always be 

1/2 

normalized [16J, through multiplication by (\tp), \ip)) ■ 

The second and third general properties of the scalar product imply that 
{a\ip), |0)) = afi (l^), |0)). This means that the eigenvalues of a self-adjoint 
bicomplex operator (associated with an eigenket that is not in the null cone) 
are hyperbolic numbers [H]. 

The simplest possible form of a bicomplex extension of the canonical 
commutation relations seems to be imbedded in (d). This assumption entails 
that none of the operators Xi and Pk belongs to the null cone. Assumption (e) 
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implies that H is not in JVC. Indeed if ifi vanished, for instance, H could 
not have an eigenvalue with £j ^ 0. Making use of the self-adjointness of 
Xi and Pk and the properties of the scalar product, one can show that the 
bicomplex number £ introduced in (d) is in fact a hyperbolic number |15j . 

Assumption (f ) is not really necessary. In effect it restricts the eigenvalues 
we will obtain to the discrete spectrum. Assumption (g) can be seen as 
contributing to the specification of the scalar product. That assumption is 
not necessary to derive eigenvalues through the algebraic method, but it is 
needed to give structure to the T-module of eigenkets. In the differential 
equation method that will be used in sect. [5J (g) will in fact be derived. 

We will now show that, without loss of generality, £ can be taken in O+Q 
To do this, we show that a simple rescaling of Xi and Pk transforms the 
problem specified in assumptions (a) to (g) into an equivalent one, but with 
£ in D+. 

For s — 1, 2, let aj and /3-$ be in C(i 1 ) and nonzero. Define X[ and P' k so 
that 



We have 



Xi := (aie 1 + a 5 e 2 ) X[ and P k : = (/^ + /%e 2 ) P' k . (15) 

(|0), X^)) = {a T e 1 + a 2 e 2 }(|0>,X^>), 
|V)) = {afe, + cF 2 e 2 } (X^), |^» . 

Since Xi is self-adjoint, the left-hand sides of these two expressions are equal. 
The right-hand sides must also be equal. Using that and requiring X[ to be 
self-adjoint, we find 

{aje 1 + Q5e a }(|0),X t qV)) = {afe 1 + a^e 2 } (](/)), X^\ip)) . 

Since X- is not in JVC, one can always find kets \<p) and \ip) such that 
(|0), X-lip)) is not in AfC. Therefore = of. By a similar argument, 

Ps=Js- 



Let us now substitute (15) into (11). In the idempotent basis we have 

-1/2' 



H 



/ 2 



(16) 



2 D + differs from D + introduced in assumption (c) in that in 
ai and «2 are excluded. 



vanishing values of 
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Sums over indices like i, j and k run from 1 to 3, whereas the idempotent- 
basis index s runs from 1 to 2. We want H to have a form similar to (11 ) or, 
again in the idempotent basis, 



-1/2* 



(17) 



Comparing (16) and (17), we see that fj! = )u//3|. This holds for both s = 1 
and 2. Hence /3| = /3| or, equivalently, /3j = ±/%. Comparing again (16) 
and (17), we have 



-1/2 



Ze 2 



^2 a ^'s 



-1/2 



Once more, the only way (Ze 2 )' can be a real positive number is if a| 
or = iasj- 

Now assumption (d) and definition (15) allow us to write 

[X,, P k ] = i^Sik + &e 2 ) I = (a T /3 Tei + a 2 /%e 2 ) [X' t , P' k ] . 
This implies that 



[XlP' k ] = i 1 h5 l , 



1 
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e 2 J =: ijMifc (^e x + £~e 2 ) I. 



Therefore, we can always choose a-g and (3$ so that £~ and £~ are real and 
positive. Moreover, we can rescale £~ to be 1, but we cannot in general rescale 
£~ so that £~ = We conclude that if can always be written as in (11), 
with the commutation relations between Xi and Pk given by 



[X h P k ] = i-JiSik (£ie x + ^e 2 ) /, fj, & e 



(18) 



3.2 Eigenvalues of H 

In this section, eigenvalues of H are obtained through Pauli's algebraic 
method. This shows the unicity of the eigenvalues obtained, under assump- 
tions (a)-(g). 
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Just as in standard quantum mechanics, we define the bicomplex angular 
momentum operator L := R x P . In terms of the Levi-Civita symbol we 
have 



Li — €ijkXjPk — < EijkXj s Pks \ 

jk s v ji ) 



(19) 



jk s K jk 

The bicomplex Runge-Lenz vector is defined as 

A := — (P x L - L x P) - Ze 2 — , 
which is equivalent to 

A i = > A - I > ' tijkPjsL ks - i^P* ) - Ze 2 ^ }> e s . (20) 



^ ] S — ( ^ tjjkPjsLks ~ l\fo£,sPis 
s [ V \ jk 



Ra 



Let us write r] := h£. We then observe that the commutator of Xi S and 
Pj S , as well as the definitions of L j S and Aj a , are the same as the ones in the 
standard quantum-mechanical case, except that K is everywhere replaced by 
r]s. By an argument identical as in the standard case [2"Tl [22] , we therefore 
obtain the following commutation relations and properties: 

[A is ,H s } = 0, (21) 

[Li s , Aj s ] = ij 7 ?? ^ tjjkAks, (22) 
k 

[Ais, Aj g ] = H s ^ ^ tjjkLks i (23) 

^ ] LjsAjg = = ^ ] Aj S Lj S , (24) 

i i 

^=(^ 2 )^F s {L 2 +?7 §}. (25) 

Here L 2 = J2iL 2 s and A 2 = X^^L- By properties of the idempotent basis, 
operators H, Li, and A,- satisfy similar relations as (21)-(25), with the index 



s deleted. Note that, because of (23), operators H s , L is , and Aj a do not 



make up a Lie algebra. They do, however, generate an infinite-dimensional 
one. 
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To avoid having to work with (the bicomplex generalization of) an infinite- 
dimensional algebra, we will restrict the action of H, Li, and Aj on the 
module Te corresponding to a given eigenvalue E of H. By assumption (f), 
J~e is finite-dimensional. Assuming as in (e) that E is not in J\fC, we define 
three operators Ai acting on T E as 



A, 



IE 



A,. 



(26) 



Equations (21)-(23) (without the index s) imply quite straightforwardly that 



A h H 

Li , Aj 

Ai , Aj 



0. 



Making use of (24) and ( |25| ), we easily obtain 

LiAi = = Y, Ai L i 



and 



(27) 
(28) 

(29) 



(30) 



(31) 



To arrive at explicit values of E, it is useful to construct six operators Fi 
and d as 



E. :-- 



1 



Li — Ai ) and d :- 



1 



Li + As). 



(32) 



Since Lj and Ai commute with H, so do F^ and Gj. The commutation 
relations of the latter are given by 



[Fi,Fj] = i^^eijkFk, 
k 

[Gii Gj] = 1-iT] €ijkGk, 
k 

[FuGjI^O. 



(33) 
(34) 
(35) 
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This means that the Fi commute with the Gj, but the Fi (and the Gj) have 
with themselves the same commutation relation as the bicomplex angular 
momentum operators. Since Fi and Gj commute, F 2 commutes with G 2 . 
With (30) and (32), one can show that F 2 — G 2 = and, therefore, that 
their eigenvalues are equal. 

If we project (33) on the idempotent basis, we find that for each s the 
Fi S have the commutation relations of standard angular momentum (with 
h replaced by rf§). The eigenvalue of F 2 is therefore equal to f$(fs+ 1) f] 2 , 
where /j is a nonnegative integer or half integer [26] . The eigenvalues of F 2 
(and of G 2 ) are therefore equal to 



(36) 



where of course / = f^e 1 + j^i- 

If we substitute d32l in pTj) we get 







2Y J F l F l + 2Y J G 3 G J+ r 1 2 + ^(Ze 2 ) 



For the eigenvalues this entails that 

= (2/ + l)V + 
or, in the idempotent basis, 



2\ 2 



IE 



(Ze 2 ) 



= W(2/, + l) 2 ,i + ^JZe 2 ) 2 



But then 



(2/,+ l)M+^(^e 2 ) 



-,2\ 2 



2£~ 



for both s — 1 and 2. This yields for E-g 

fi(Ze 2 f ^{Ze 2 f 



2^,4 



2r / |(2^+l) 2 "" 2ft 2 £f(2/ ? +l) 2 
where n^- is a positive integer. This means that we can write 

/xZ 2 e 4 r /iZ 2 e 4 



2fr 2 £fn. : 



2^,2 ' 



-En 



2fr 2 £% 2 



E 



2n 2 ^|n| 



(37) 
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This coincides with the standard Coulomb potential energy levels 



if 



and only if £j 



£2 and ni 



n. We call n = n^e 1 + n^&i the 



bicomplex principal quantum number. In the hyperbolic representation. 
//Z 2 e 4 



Ah 2 



+ 



^-2 -2 
s 2 2 



We point out a formal symmetry satisfied by (38). Since j : 



(38) 



J 2- 



we can see that a/j = e i + ii e 2- But then £v/j = ^e-,^ + i^^- This 
immediately implies that 

Re{E n ,£) = Hy^^Vl) and Me^^Vl) = . 

Note, however, that £ a/j is not a hyperbolic number. 



4 Eigenfunctions of 

In this section, we define a coordinate-basis representation for the bicomplex 
operators Xi and Pk- We then show that they, as well as H, are self-adjoint 
with respect to the scalar product defined in sect. |2.2| Finally, solving the 
hamiltonian eigenvalue equation in the coordinate basis, we recover eigenval- 



ues given in (37) and obtain the hyperbolic Coulomb potential eigenfunctions. 



4.1 Coordinate-basis representation 

We begin by constructing a representation of Xi and Pk on the space 
of bicomplex square- integrable functions on IR 3 . Letting r denote the triplet 
(xi, x 2 , x 3 ), we define the action of Xi as 

XJ(r) := xj(r) . (39) 

This implies that the Xj commute two by two. A function F(R) acts on / 



MS 



F(R)f(r) := F(r)f(r) 



(40) 



For the action of Pk we write, in a rather straightforward extension of the 
standard case 



P k f(r) := -hf*£-M = -h.l-m 



(41) 
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where £ = + ^ e 2 an d the £j are positive real numbers]^] Clearly the P k 
commute two by two. Moreover, by letting both sides act on an arbitrary 
function /, one easily shows that 

[X u P k } = i lV 5 lk I. (42) 

Let / and g be in J^. Clearly (Xif,g) = (f,Xig), so that Xi is self- 
adjoint. For P k we have 

(P k f,g)-(f,P k g) 




To obtain the last equality, we have restricted the space of functions to those 
that vanish at infinity and on these functions, P k is self-adjoint. The proof 
that (Hf,g) = (f,Hg), and therefore that H is self-adjoint, is straightfor- 
ward. 



4.2 Wave functions 

The bicomplex quantum Coulomb potential problem consists in solving the 
three-dimensional eigenvalue equation 



Htl) E {r) = Ei/} E (r) 



(43) 



for if given by (11). Making use of (40) and (41) we can write more explicitly 

Ze 2 



z/i r 



i[> E (r) = Eijj E (r) . 



(44) 



We now write £, E, and ip E in the idempotent basis. Equation (44) 
becomes 



E 



2/i 



Ze 2 
r 



ih V 2 + — + Es\(ij E Ur)e s = 0. 



(45) 



3 Strictly speaking, P k should be denned on a subset of J-3, made up of suitably diffcr- 
entiable functions. We'll come back to this in sect. [6] 
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Clearly, each coefficient of e s must separately vanish. Writing the laplacian 
in spherical coordinates and making use of the standard expression of the 
angular momentum operator we get 



2/x 



1 d 



7" ^ ()lT V C^T 1 







l_Lf 

is 



Ze 2 

+ — + (^)j = 0, (46) 



for s — 1 and 2. Now we know from standard quantum mechanics that the 
spherical harmonics are eigenfunctions of the square of the angular momen- 
tum, that is, 

L 2 t Y l?mt = k(k+l)4Y hmt , 
with —Is < rris < h- Looking for solutions of the form 



(47) 



(48) 



and using (47), we get for s = 1,2 

~k(k+l) 2/2 fZe 2 



1_ d_ 

r 2 ^ r 



r -rUs(r) 
dr 



4 



u^{r) = 0. (49) 



Just as in the standard case, the suitably normalized solutions of (49) are 
given bv [M 125] 



f 2Z \ 3 (ng-fg-i; 
V^sd/ 2n s -[(n 3 -+ fe) 



1/2 



(50) 



where < nj, the L^£^i_! are Laguerre polynomials and 



Cs ■= r, 

n-sdos 



4_ 

- lie 2 ' 



(51) 



Each solution corresponds to a value of E-g given by —fiZ 2 e 4 /2r] 2 n 2 . Thus 
we recover the eigenvalues (37). The bicomplex wave functions in (43) can 
be written as 



ipnim(r) = u nl {r)Y lm (6 , (p) = ^u nth (r)Y h 



(52) 
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For £j and £5 fixed, any sextuplet (n-j, n^, Z-j, ^f, defines an eigen- 
function of H. All functions with the same (n-j, n-g) correspond to the same 
eigenvalue. A general eigenf unction of H can therefore be written as 

ns— 1 

EE E ^W^WM**. with Q«6T. (53) 

s £j=0 m-g=—l-g 

Along the way, we have introduced the bicomplex orbital quantum number 
I = l^e 1 + l^e 2 and the bicomplex magnetic quantum number m = m^e 1 + 
m^e 2 . When = n^, the number n is real, and similarly with / and m. 



5 Graphical representation 



It is interesting to plot some of the eigenfunctions (52) explicitly. We will 
concentrate on the radial part. (Note that if Z-j 7^ Z5, the spherical harmonics 
will never factor out as real functions of 9 and 0.) We write 

u n i{r) := it„ T i T e 1 + u n ^e 2 . (54) 

It is particularly instructive to use the decomposition of {e 1; e 2 } in terms of 



j given in sect. 2.1 We can then rewrite (54) as 



UnlVT 



s 

+ \ E (-!) S+1 V^^^SUiC^j. (55) 



where \ ul , fr d is the normalization constant. 

y rl s i s s 

It will be useful to define, as in ([3]), 

Ct = x c + y c , (2 = x (- Vc, ( 56 ) 

and consider the three functions M.e(u n i), Hy(w n z) and \u n i\ as depending on r 
through the two variables and y^. With £j fixed, graphical representations 
of these functions can be obtained by assigning specific values to rtg and 
Figure [l] shows the real part, the hyperbolic part and the real norm of ( 55 ) 



in the specific eigenstate corresponding to = 8, = 7, = 6, % = 5. 
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0.035 




150 200 



(a) ao /2 Me(u 8 7 65 ). (b) ao /2 Hy(u 8765 ). (c) a^ 2 |u 8765 |. 

Figure 1: Coulomb potential radial function in a specific eigenstate, for £5 — 
0.50 (red), 0.62 (green), 0.74 (blue), 0.86 (yellow) and 0.98 (purple). 

Here £j = 1 and we have shown the curves obtained for different values of £5, 
from 0.50 to 0.98 by intervals of 0.12. The independent variable p is defined 
as 2Zr/ao and the function plotted is indicated under each graph. 

In fig. |2j we go beyond the representation of eigenf unctions of H and 
extend and to genuinely independent variables. For £ fixed, this is 
equivalent to considering r in the hyperbolic plane. This allows for interest- 
ing surfaces to emerge, connected to a new class of polynomials in two real 
variables. We take £1 = 1 = £2 an< ^ ^ n » = 25 and Zg- = 12 for s = 1, 2. In 
these plots, a cut at = makes the hyperbolic part of u n i vanish. The real 
part and real norm then coincide with the ones in the standard case. 

There is another equivalent way to generate the radial surfaces in the 




(a) ao /2 Me(M 2525 i2i2). (b) ao /2 Hy( M 2525i2i 2 ). (<0 4 \u 25251212 \ 2 . 

Figure 2: u n i as a function of two independent variables. 
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hyperbolic basis {1, j} in the standard case. We can write eigenfunctions (54) 
as 

where £ := x^ + y^j, ( := x^ + y^j and 

MO :=C'^-i(0 =Mtm(x ( ,y<;))+My(e nl (x y c ))j, 

where 

^i(C):=E L K-i^)e.. 



Using the correspondence between the idempotent and hyperbolic bases, one 
can easily show that 

r 3 = \{(^+ %r 3 + - %r 3 } + ^ {(^ + %)~ 3 - - %)~ 3 } j> 

which we define as + y^j . It is not difficult to show that the exponential 
transforms as 



-C/2 



e-^ 2 {cosh|-jsinh|} 



This suggests that we can explicitly write 



4 cosh y ~ y t 



t sinh 



le(£ n i(x ( ,y c )) 



+ (y'j: cosh y - 4 sinh ^) ^M^O fc)) 



Hy(u 



4 cosh ^ - 2/£ sinh 7T ) ^(^ni^C' fc)) 



+ ^ cosh y - 4 sinh f") Re(4j(a: c , %)) 



With a symbolic computation software, we can generate the polynomial £ n i{() 
for any positive integers and l-g < rig — 1, for s — 1,2. We then take 
n — 25, I — 12, = 1, y% = (implying that x'^ = 1 and y'^ — 0), compute 
4ii(C)> transform £ into x^ + y^j and separate the real and hyperbolic parts 
to get M.e(£ n i(x^,y^)) and M.y(£ nl (x^,y^)) explicitly. If we plot the associated 
a^ 2 M.e(u), Oq Hy(tt) and \u\ 2 , we recover the results shown in fig. [2 



18 



6 Discussion 



We have solved the eigenvalue equation (43) for the discrete spectrum of the 



Coulomb potential hamiltonian (11 ) in the framework of bicomplex numbers. 



The continuous spectrum could also be worked out along similar lines. The 



eigenvalues corresponding to the discrete spectrum are given in (37) and the 
eigenfunctions in (52). Note that if £j- = 1 = the standard wave functions 
can be recovered by letting rv^ = = and = m^. 

It is instructive to investigate the orthogonality properties of the eigen- 
functions (52). Making use of definition ([9j of the scalar product, one can 
write 



[ipnlm, tpn'l'm' 



i ) nlm( r )' l l } n'l'm'(r)dr 



It is well-known [2lJ EH] that the standard Coulomb problem eigenfunctions 
are orthonormal in all indices. This implies that 



(58) 



From (58) we can draw two conclusions: 



1. The eigenfunction ip n i m is normalized. Indeed 



(59) 



2. If E n — E n i is not in the null cone, then ipnim and Vvz'm' are orthogonal. 
Indeed from (37) we see that E n — E n > is not in the null cone if and 
only if 7^ n~ for s = 1,2. But then 5 n ~ n i 



for 



1, 2 and the 



orthogonality follows from (58). 



Let us consider the set of all finite linear combinations of eigenfunctions 
ipnim, with bicomplex coefficients C n i m . It is straightforward to show that 
this set makes up a T- module, which we denote as AA. Defining Xi and Pk 



as in (39) and (41), one sees that these operators are well-defined on Ai. 
Moreover, it is not difficult to show that properties (a)-(g) in sect. 3.1 are 
satisfied in A4, thus proving their consistency. 
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Of course, M. is not a Hilbert space, since the restriction to finite lin- 
ear combinations entails that it is not complete. It is well-known that the 
Coulomb potential eigenfunctions in standard quantum mechanics, i.e. the 
V'njJjms; make up an orthonormal system in the Hilbert space L 2 (IR 3 ) [5]. 
From [T7j, one concludes that the il> n i m make up an orthonormal system in a 
bicomplex Hilbert space L 2 (R 3 ) © L 2 (IR 3 ). 

We close this section with a result that we prove with the notation of 
the Coulomb problem, but that clearly holds more generally. Let U be a 
bicomplex linear operator, acting on Ai, that commutes with H. Then 

{E n - E n i) (ll>Tdm,Ull>n'l'm') = 0. (60) 

The proof is straightforward: 

= {ipnlm, [H, U] Ipn'i'm') 

= {ipnlm, {HU -UH} tpn'l'm') 
= {ipnlm, HUtjj n >ii m i) — {vfj n i m ,UHvjj n iii m i) 
= {U* H*1p n i m ,1p n /i/ m i) — {lpnlm-,UHlp Tl iii m i) 
= E n {U*1p n l m , Ipn'l'm') ~ E n i {lp n l m , Ulp n 'l' m ') 
= {E n — E n i) {ipnlm, U Ipn'l'm') ■ 

This means that if E n — E n i is not in the null cone, then {ip n i m , Utp n 'i' m ') van- 
ishes. In other words, Uip n 'ir m > is a linear combination of functions associated 
with eigenvalue E n >. 

7 Conclusion 

We have shown that, just like the quantum harmonic oscillator problem [15] . 
the quantum Coulomb problem can be solved in the framework of bicom- 
plex numbers. We have obtained the eigenvalues and eigenfunctions of the 
bicomplex hamiltonian, and have shown that the eigenfunctions make up an 
orthonormal system in a bicomplex Hilbert space. The question is still open 
whether the constants £j and £g can be given a physical interpretation. In 
any case, it is likely that the mathematical properties of the functions we 
introduced can be fruitfully studied for their own sake. 
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